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ABSTRACT 


This  work  simulates  an  acousto-optic  signal  processing  technique 
that  can  be  used  to  filter  radio-frequency  (RF)  electronic  signals.  Of 
particular  interest  is  the  removal  or  excision  of  narrowband 
interference  from  broadband  signals.  The  program  developed 
simulates  the  performance  of  a  binary  phase  shift  keying  (BPSK) 
signal  with  a  narrowband  interference.  The  width  and  depth  of 
notches  that  model  an  excisor-detector  array  are  variable.  In 
addition  there  is  a  capability  to  move  the  notches  in  order  to  place  a 
notch  over  the  interference  for  maximum  filtering  capability.  A 
summary  of  several  possible  excisor  techniques  is  presented. 
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DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in 
this  research  may  not  have  been  exercised  for  all  cases  of  interest. 
While  every  effort  has  been  made,  within  the  time  available,  to 
ensure  that  the  programs  are  free  of  computational  and  logic  errors, 
they  cannot  be  considered  validated.  Any  application  of  these 
program  without  additional  verification  is  at  the  risk  of  the  user. 
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I .  INTRODLCTION 


A-  BACKGROUND 

Many  companies  are  designing,  developing,  and  field-testing 
Acousto-optics  Spectrum  Analyzers  (AOS A)  for  communications, 
radar  and  electronic  warfare  surveillance.  In  many  such 
applications,  acousto-optics  is  not  merely  a  good  choice —  it  is  the 
only  choice.  [Ref.  1] 


Figure  1  -  Acousto-optic  spectrum  analyzer 


As  an  introduction  to  acousto-optic  spectrum  analyzers  consider 
the  system  in  Figure  1.  Coherent  light  illuminates  a  one  dimensional 


spatial  light  modulator,  in  this  case  a  Bragg  cell.  The  light  diffracted 
by  the  cell  is  collected  by  the  lens  and  focused  onto  a  photodiode 
array.  The  light  amplitude  at  the  photodetector  is  the  spectrum  of 
the  information  in  the  Bragg  cell.  This  serves  to  illustrate  the 
simplicity  of  one  optical  technique  for  calculating  spectra.  [Ref.  2] 

The  problem  of  rejecting  relatively  narrowband  interference 
signals  is  becoming  more  important  as  the  bandwidth  of  digital 
communications  systems  increases.  Acousto-optics  is  showing 
promise  in  this  area  [Ref.  1].  An  optical  implementation  of  this 
technique  is  shown  in  Figure  2. 

In  order  to  implement  an  optical  excisor  we  need  more  elements 
interacting  in  the  AOSA  basic  architecture  than  was  presented  in 
Figure  I.  As  showed  in  Figure  2,  the  signal  s(t)  is  applied  to  a  Bragg 
cell  (1).  The  cell  is  illuminated  with  coherent  light  (2)  and  the 
spectrum  of  the  signal  is  formed  on  the  photodetector  array  (3),  The 
light  at  each  photodetector  has  an  amplitude  proportional  to  the 
amplitude  of  the  corresponding  acoustic  frequency  in  the  Bragg  cell. 
A  given  acoustic  frequency  component  acts  as  a  moving  diffraction 
grating;  for  this  reason  the  light  diffracted  by  it  is  Doppler 
frequency-shifted  by  the  frequency  of  the  acoustic  signal.  The 
acoustic  signal  component  can  be  reconstructed  by  down-converting 
the  diffracted  light  to  the  acoustic  frequency.  This  is  achieved  by 
adding  the  optical  signal  to  a  reference  beam  (4).  The  resulting 
combined  beams  fall  on  a  photodetector  array.  The  photodetector 
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produces  an  electrical  signal  with  a  frequency  equal  to  the  difference 
between  the  frequencies  of  the  two  optical  inputs.  This  results  in  the 
recovery  of  the  acoustic  frequency  components  and  hence  the 
recovery  of  the  input  radio  frequency  component. 

The  optical  excisor  continuously  forms  the  spectrum  of  a 
broadband  electrical  signal  and  reconstructs  an  exact  replica  of  the 
input  signal  from  that  spectrum.  The  important  feature  of  this 
process  is  that  the  spectrum  is  spatially  distributed  and  is  physically 
available  at  the  detector  plane  for  removing  an  undesired  portion  of 
the  spectrum.  For  example,  the  excision  of  narrowband  interference 
from  broadband  signals  is  possible. 


Figure  2  -  Optical  excisor  architecture 


Two  significant  characteristics  of  an  optical  excisor  are:  1)  it  has 
an  RF  signal  input  and  an  RF  signal  output,  2)  it  is  a  predetection 
signal  processor  producing  a  continuous  coherent  output  [Ref.  3]. 
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Figure  3  -  Optical  excisor  characteristic 

B  -  OPTICAL  EXCISION 

If  we  have  two  input  signals  that  are  sufficiently  different  in 
frequency  to  be  resolved  in  the  back  focal  plane,  a  spatial  light 
modulator  can  be  placed  in  this  location  and  used  to  turn  off  either 
one  or  both  of  the  input  signals  [Ref.  3]. 

Another  technique  would  be  to  block  the  light  of  one  signal  by 
using  an  aperture  or  small  wire  in  the  focal  plane,  thereby  allowing 
passage  of  only  the  desired  signal.  [Ref.  3] 

A  broadband  signal  that  has  a  narrowband  interference  signal  in 
its  spectrum  can  be  applied  to  the  input  of  the  optical  excisor.  The 
narrowband  signal  will  form  a  small  spot  of  light,  while  the 
broadband  signal  will  extend  over  a  relatively  large  area.  Using  a 
small  wire,  the  strength  of  the  interference  can  be  significantly 
reduced  while  only  blocking  a  small  portion  of  the  broadband  signal. 
The  bigger  the  wire,  the  more  the  interference  will  be  blocked,  and 
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therefore  the  more  it  will  be  attenuated  at  the  output.  However,  a 
large  wire  will  also  remove  more  of  the  broadband  output  signal  and 
perhaps  begin  to  cause  significant  degradation  in  the  broadband 
output  waveform.  The  task  for  the  designer  of  an  optical  excisor  is  to 
optimize  the  transfer  function  to  achieve  a  useful  amount  of 
attenuation  of  the  narrowband  interference  without  removing 
significant  portions  of  the  broadband  signal.  [Ref.  3] 

An  alternative  excision  approach  is  to  place  a  material  such  as 
photochromic  glass  or  photodichronic  crystal  in  the  focal  plane. 
These  materials  will  darken  in  an  amount  proportional  to  the 
irradiance  of  the  incident  light  (much  like  photographic  film,  but  in 
real  time).  Narrowband  interference  signals,  because  of  their 
relatively  high  intensity,  will  be  attenuated  automatically  [Ref,  3]. 

As  can  be  seen,  using  different  excision  methods  in  the  optical 
Fourier  transform  plane  will  create  excision  notches  with  different 
properties. 

The  following  is  a  summary  of  the  different  experiments  found 
in  references  which  illustrate  various  techniques  for  excision  as  well 
as  their  results. 

Thin  opaque  wires  of  known  diameters  placed  in  the  transform 
plane  can  create  narrow  excision  notches.  Using  thin  bare  copper 
wire  in  the  transform  plane,  it  blocks  a  portion  of  the  diffracted 
signal  spectrum,  thereby  creating  an  excision  notch.  In  different 
experiments  it  was  found  that  the  attenuation  achieved  was  about  40 
dB  with  inexpensive  lenses.  [Ref.  6] 
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A  Pockel's  Readout  Optical  Modulator  (PROM)  spatial  light 
modulator  can  be  used  in  a  fashion  which  automatically  attenuates 
narrowband  interferences.  The  PROM  has  nonlinear  clipping 
characteristics  which  attenuates  intense  signals,  while  weaker  signals 
pass  unaffected.  In  fact,  a  strong  signal  can  be  attenuated  below  the 
level  of  a  weak  signal.  Broadband  signals  should  be  relatively 
unattenuated  compared  to  more  powerful  narrowband  signals 
because  their  power  spectrum  distribution  in  the  optical  Fourier 
transform  plane  will  be  below  the  clipping  threshold  of  the  PROM. 

The  PROM  optical  clipper  is  intended  to  be  utilized  as  an  image 
storage  device  which  is  written  with  blue  light  and  read  out  with  red 
light.  The  experiment  reported  in  reference  3  utilized  the  red  HeNe 
laser  light  for  both  writing  and  reading  functions.  A  strong 
narrowband  interference  signal  will  create  a  bright  red  spot  in  the 
transform  plane  which  effectively  acts  to  write  a  dark  spot  on  the 

PROM.  This  dark  spot  prevents  the  narrowband  interference  from 
being  passed  to  the  photodetector.  Thus  the  strong  narrowband 
interference  signal  effectively  attenuates  itself.  [Ref.  3] 

Reference  3  showed  that  the  PROM-based  excisor  has  attenuated 
a  narrowband  signal  by  28  dB  relative  to  the  broadband  pulsed 

signal.  It  demonstrated  encouraging  results,  but  the  utility  in  a 
practical  system  is  limited  by  the  fact  that  a  PROM  is  an  image 
storage  device  which  must  be  cycled  through  an  erase  phase 

followed  by  a  read/write  phase. 
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Variations  of  the  PROM  are  being  considered  which  will  avoid 
the  image  storage  properties  of  the  current  PROM  so  that  time¬ 
framing  and  cycling  problems  may  be  avoided.  [Ref.  3] 

A  Programmable  lead  lanthanum  zirconate  (PLZT)  ceramic 
modulator  can  also  be  used  in  the  transform  plane.  A  programmable 
PLZT  ceramic  light  modulator  will  allow  the  filter  magnitude 
response  to  be  controlled  by  electronic  means.  By  applying  a 
transverse  electric  field  to  the  PLZT  ceramic,  the  polarization  of  the 
transmitted  optical  beam  can  be  rotated.  By  placing  the  PLZT 
ceramic  between  crossed  optical  polarizers,  the  intensity  of  the 
transmitted  optical  beam  can  be  electronically  controlled.  Spatial 
control  of  the  transmitted  optical  intensity  can  be  obtained  by  having 
separately  addressable  electrodes  across  the  aperture  of  the  PLZT 
ceramic.  [Ref.  3] 

The  results  of  the  PLZT-based  optical  excisor  showed  that  18  dB 
of  excision  could  be  obtained  while  an  opaque  wire  having  a  similar 
notch  width  obtained  23  dB  of  excision  depth.  [Ref.  3] 

Another  form  of  excision  is  due  to  the  spaces  between  detectors 
in  a  detector  array  to  perform  the  optical  excision.  In  the  frequency 
domain,  they  are  represented  by  a  set  of  notches  that  are 
periodically  spaced.  The  width  of  the  notches  is  proportional  to  the 
spaces  between  the  detector  elements.  The  depth  depends  on  the 
width  and  the  diffraction  effects  (A  similar  model  also  applies  when 
a  linear  array  of  fiber  optic  cables  leading  to  individual  detectors  is 
used  to  measured  the  light  properties).  It  is  this  latter  form  of 
excision  that  we  wanted  to  study. 
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C-  THESIS  OBJECTIVE 


The  objective  of  the  research  described  in  this  thesis  was  to 
develop  a  computer  program  which,  through  simulation,  would 
provide  a  capability  to  evaluate  the  filtering  ability  of  gaps  between 
detectors  in  a  detector  array  versus  narrowband  interference  in  a 
broadband  signal. 

Using  the  Naval  Postgraduate  School's  library  and  computer 
facilities,  we  wanted  to  develop  a  computer  program  based  in 
MATLAB386  [Ref.  4]  using  a  Compaq  386/20  computer. 

The  specific  goals  were  to  evaluate  the  distortion  on  the  original 
signal  due  to  the  filtering  process  and  to  examine  the  notch 
dimensions  (width  and  depth)  for  their  effect  on  the  filtering 
capability. 

Chapter  2  serves  as  a  brief  description  of  the  optical  excisor 
system.  This  chapter  describes  an  acousto-optic  signal  processing 
technique  that  can  be  used  to  filter  radio-frequency  electronic 
signals.  A  discussion  of  the  mathematical  model  chosen  for  the 
simulation  is  given  in  Chapter  3.  The  detailed  description  of  the 
computer  program  is  presented  in  Chapter  4,  showing  the 
information  available  during  the  simulation.  Chapter  5  contains  the 
simulation  results  and  an  evaluation  of  the  program,  and  finally 
Chapter  6  contains  concluding  remarks  and  suggests  areas  of  future 
research. 
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II.  GENERAL  DESCRIPTION  OF  THE  SYSTEM 


The  optical  excisor  separates  the  input  signal  s(t)  into  each  of  its 
Fourier  components,  processes  the  magnitude  and  phase  of  each 
Fourier  component  and  finally  sums  the  modified  Fourier 
components  into  a  filtered  signal.  A  general  description  of  an  optical 
excisor  can  be  based  on  the  architecture  shown  in  Figure  5. 

The  electrical  input  signal  s(t)  typically  is  derived  from  a 
broadband  antenna,  downconverter,  and  amplifiers.  This  signal  is 
applied  to  a  Bragg  cell  illuminated  by  a  coherent  light,  illustrated  in 
Figure  4. 


Figure  4  -  Bragg  diffraction  in  a  spectrum  analyzer 
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The  acousto-optic  modulator  is  a  block  of  transparent  material 
with  a  piezoelectric  transducer  bonded  to  one  end.  The  other  end  is 
tapered  to  trap  the  acoustic  waves  and  loaded  to  provide  an  effective 
acoustic  termination.  A  piezoelectric  transducer  is  sandwiched 
between  electrodes  and  vibrates  in  synchronism  with  the  applied 
voltage.  When  driven  by  an  electrical  signal,  the  transducer  excites 
an  acoustic  which  is  a  replica  of  that  signal.  This  acoustic  wave 

propagates  down  the  long  axis  of  the  Bragg  cell  and  produces 

compressions  and  rarefractions  in  the  material  [Ref.  2].  The 
interaction  between  coherent  light  and  the  acoustic  wave  results 
from  an  increase  in  the  index  of  refraction  of  the  material  at  a 

compression  and  a  decrease  in  the  index  at  a  rarefraction.  This 
results  in  a  modulation  in  phase  of  the  light  passing  through  the 
acousto-optic  modulator.  The  resulting  spatial  phase  distribution  can 
be  considered  equivalent  to  a  superposition  of  plane  waves  of 
various  phases  and  amplitudes.  By  using  a  lens  to  focus  the  various 
diffracted  plane  waves,  each  wave  can  be  separated  and  the 

amplitude  and  phase  of  each  spectral  component  of  the  original  input 
signal  can  be  measured.  So  far  the  description  is  similar  to  an  AOSA 
[Ref.  31. 

In  addition  to  diffracting  the  plane  waves  based  on  the  spectrum 
of  the  input  signal,  the  acousto-optic  modulator  also  Doppler  shifts 
the  diffracted  plane  waves.  The  amount  of  Doppler  shift  is  exactly 
equal  to  the  frequency  of  the  input  signal,  and  therefore,  if  the 


amount  of  Doppler  shift  is  measured,  the  frequency  of  the  input  can 
be  determined.  [Ref.  3] 

An  easy  way  for  measuring  an  optical  Doppler  shift  is  to  add 
some  unshifted  light  from  the  optical  source  and  process  the  summed 
signal  electronically.  The  resulting  beat  signal  frequency  will  be 
equal  to  the  frequency  difference  between  the  two  signals  and  hence 
be  equal  to  the  amount  of  Doppler  shift.  [Ref,  3] 

The  unshifted  reference  beam,  in  Figure  5,  is  derived  from  the 
spatially  coherent  optical  source  using  optical  splitters  (partially 
reflecting  mirror).  The  light  reflected  by  this  beam  splitter  is  set  to 
mirror  1  which  then  reflects  the  light  to  the  acousto-optic  modulator. 
The  light  transmitted  by  the  beam  splitter  is  reflected  by  mirror  2 
up  to  a  second  beam  splitter,  reflecting  some  of  the  reference  beam 
to  the  photodetector  array.  The  second  beam  splitter  also  passes 
some  of  the  light  diffracted  by  the  acousto-optic  modulator.  By 
adjusting  the  mirrors  and  beam  splitters,  both  the  reference  beam 
and  the  diffracted  beam  appear  to  come  from  the  same  direction  and 
are  added  together  by  the  beam  splitter  before  the  detector  array. 
By  using  a  spatial  light  modulator  in  the  transform  plane  and 
selectively  passing  only  the  desired  portions  of  the  spectrum, 
interference  signals  can  be  rejected  [Ref.  3].  The  photomixing  and 
signal  summing  by  the  detector  reconstructs  the  time  domain 
waveform  of  the  desired  portion  of  the  input  signal  and  provides  a 
coherent  continuous  output  i(t). 


The  five  main  components  of  the  system  presented  are: 

(1)  spatial  coherent  light  source; 

(2)  one-dimensional,  acousto-optic  modulator; 

(3)  lens  for  Fourier  transformation; 

(4)  a  frequency  plane  spatial  modulator,  and 

(5)  the  photodetector  array. 

The  present  work  considers  the  detector  array  acting  as  a  spatial 
light  modulator  in  the  frequency  plane;  for  this  reason,  from  now  on, 
it  will  be  called  as  "excisor-detector"  array.  This  technique  is  of 
particular  importance  in  rejecting  narrowband  interference  in 
wideband  communications  systems.  By  using  the  spaces  between 
detectors  as  notches  in  the  frequency  domain,  one  can  control  the 
frequency  rejection  to  match  the  signal  environment.  The  signal 
environment  may  consist  of  a  wideband  signal  with  a  narrowband 
interference  signal  located  at  some  arbitrary  frequency  within  the 
bandwidth  of  the  desired  wideband  signal.  The  object  is  to  block  the 
narrowband  interference  signal  with  one  of  the  narrow  bandwidth 
frequency  rejection  notches  of  the  excisor-detector  array. 


III. 


THE  MATHEMATICAL  MODEL 


A  more  precise  description  of  the  operation  in  the  optical  excisor 
can  now  be  given.  The  broadband  signal  including  a  narrowband 
interfering  signal  is  s(t).  The  collimated  coherent  input  light  is 
diffracted  through  an  angle  proportional  to  the  acoustic  frequency 
present  in  the  Bragg  cell.  This  acoustic  frequency  results  from  the 
applied  electrical  signal  s(t),  where  s(t)  is  real.  The  diffracted  optical 
signal  amplitude  at  the  exit  port  of  the  cell  is  proportional  to  s(t-x/v), 
where  x  is  the  distance  from  the  input  transducer  and  v  the  acoustic 
velocity.  The  acousto-optic  (AO)  cell  has  a  finite  aperture  equal  to  T, 
where  T  is  the  transit  time  of  the  acoustic  wave  in  the  AO  cell.  This 
finite  aperture  corresponds  to  an  apodization  of  the  acoustic  signal 
such  that  the  signal  may  now  be  represented  by  s(t-x/v-T/2)  [Ref.  2]. 
Also  for  this  analysis,  the  AO  ceil  is  assumed  to  operate  over  an 
octave  RF  bandwidth  (or  less)  and  the  positive  first  order  diffraction 
is  used  [Refs,  2,  3  and  5]. 

The  optical  system  takes  the  Fourier  transform  with  respect  to 
the  spatial  variable  x. 


S(p,t) 


dx 


(1) 


where  L  is  equal  to  the  length  of  the  AO  cell  and  p  is  the  spatial 


frequency  variable.  Associated  with  each  spatial  frequency  p  is  a 
temporal  frequency  co  =  pv. 

A  change  of  variables  (where  u=t-x/v-T/2)  yields: 

c-  \  “  j(t  ~  y)p'*'  f  ^  ^  .  X  jpv  u  . 

S(p,  t)=ve  2  J  s(u)e  du 

‘  (2) 

S(pv )  =  I  s(u)  du 

‘  (3) 

Equation  (2)  is  seen  to  be  the  Fourier  transform  of  the  portion  of 
s(t)  within  the  aperture  of  the  Bragg  cell.  The  exponential  phase 
term  indicates  that  the  light  diffracted  to  position  p  is  frequency 
shifted  by  an  amount  equal  to  the  temporal  frequency  co. 

The  detectors  respond  to  light  intensity  rather  than  light 
amplitude,  for  this  reason  it  is  necessary  to  bring  in  a  reference  wave 
intersecting  the  frequency  plane.  The  reference  wave  has  a  flat 
wavefront  which  may  be  parallel  with  the  frequency  plane  or 
intersect  the  frequency  plane  at  a  small  angle  0.  This  reference  wave 
with  tilt  angle  0  may  be  written  as  Ae^'^P,  where  A  is  the  amplitude 
of  the  reference  wave  [Ref.  5]. 

The  total  amplitude  in  the  output  plane,  T(p,t),  is  the  sum  of  the 
reference  wave  and  the  transform  of  the  signal  multiplied  by  the 
filter  signal,  H(p,t),  (representing  the  spaces  between  detectors). 

T(p,t)  =[AeJ‘^P  +  v  e"-'P''^^  "  T^S(pv)]  x  H(p,  t) 


(4) 


The  corresponding  intensity  distribution  is: 


I„(p,.)  =  |Tl(p,t)| 


2 


(5) 


|T(p,t)f  =  H(p,t)^A^  +  V ^  S(pv)^H(p,t)^ 

+  2  Re[S  (p  v)^H  (p, t)^  e  j -  f 


(6) 


where  Re(x)  means  the  real  part  of  (x).  The  output  from  the  excisor- 
detector  array  is  the  summation  of  all  photodiodes;  this  corresponds 
to  an  integration  over  the  spatial  variable  p. 


/  I  o(p,  t)  dp  =  /  |T(p,  Op  dp 

=  jjH(p,  t)v  e"j''^^  “f^S(pv)  +2Aej'^PHtp,  t)|^  dp 
=  |v^  S(pv)^  H(p,  t)  dp+A^/H(p,t)^  dp 

T 

+2Av  RejH(p,  t)S(pv)  ■  2^^  dp 


(7) 


(8) 


As  stated  before,  the  s(t)  signal  is  assumed  to  be  real  and  to 
occupy  an  octave  or  less  of  bandwidth.  Under  these  assumptions  the 
first  two  terms  of  equation  (8)  do  not  contribute  in-band  frequency 
components  to  the  reconstruction  of  the  original  input  signal  [Refs.  2 
and  5]. 
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The  last  term  of  equation  (8)  represents  an  inverse  Fourier 
transform  of  the  original  Fourier  transformed  signal,  which  yields  the 
original  signal  with  time  delay  equal  to  half  the  time  aperture  of  the 
AO  cell  plus  the  delay  equal  to  the  reference  beam. 

T 

2AvRej H(p,t)S(pv)eJ>^v‘'^  "T**  dp  = 

=  2Ah(t)s(t-i  +  i) 

^  /C\\ 


Equation  (9)  represents  the  filtered  reconstruction  of  the  original 
signal. 

The  mathematical  model  described  can  be  summarized  in  Figure 
6  for  the  program  calculations.  The  test  BPSK  signal  Vi(t)  is  injected 
in  the  input.  The  quantity  Vi(f)  represents  the  FFT  of  the  input 
signal,  H(f)  is  the  transier  functic:'  of  the  detector-excisor  array,  and 
Vo(f)  represents  the  filtered  signal.  Recovering  the  filtered  signal 
back  to  the  time  domain,  the  IFFT  is  performed  to  obtain  Vir(t) 
which  is  electronically  processed  to  detect  the  original  binary 
information. 
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IV.  COMPUTER  PROGRAM  DESCRIPTION 


A  -  INTRODUCTION 

The  purpose  of  the  simulation  program  is  to  model  and  analyze 
the  effects  of  an  optical  excisor.  Of  particular  interest  is  the  removal 
of  narrowband  interference  from  a  BPSK  wommunication  signal.  This 
interference  is  removed  using  the  notches  between  detectors,  in  a 
detector  array,  as  a  filter. 

Optical  excisor  performance  is  measured  in  terms  of  how  well 
the  offending  signal  frequency  components  are  removed  from  the 
output,  and  how  much  of  the  original  signal  of  interest  is  corrupted 
through  loss  of  frequency  components  during  the  filtering  process. 
[Ref.  6] 

The  excisor-detector  array  contains  gaps  between  detectors  that 
appear  as  notches  in  the  frequency  domain  (filter).  For  our  purpose 
the  width  and  depth  of  the  notches  can  be  varied  to  simulate  optical 
and  mechanical  adjustments  and  the  effects  produced  on  the  test 
signal  can  be  observed. 

The  modeling  is  intended  to  be  simple  enough  as  not  to  obscure 
the  fundamental  ideas  but  complete  enough  to  include  important 
effects  and  avoid  incorrect  conclusions. 


B-  PROGRAM  STRUCTURE 


The  program  was  divided  into  eight  different  blocks  that 
represent  various  stages  in  the  process.  Each  of  these  blocks 
performs  specific  functions  presented  in  Table  I. 

A  simplified  diagram  of  the  program  is  presented  in  Figure  7 
which  gives  the  general  structure  of  the  system. 


TABLE  I  -  BLOCK  FUNCTIONS 


BLOCKS 


FUNCTIONS 


•  INP_B.ni 

•  M.BPSK.m 


•  FILTER.m 

•  S.BPSK.m 


•  FFT_B.m 


•  S_FIL.in 


•  IFFT  B.m 


1-  Lets  you  input  the  data  for  the  simulation. 

1-  Master  Block  executes  other  sections. 

2-  Controls  the  loops  to  simulate  the  movement  of 
the  Excisor-detector  array. 

1-  Builds  the  filter. 

2-  Generates  the  frequency  axis. 

1-  Creates  the  BPSK  signal. 

2-  Creates  the  interference. 

3-  Generates  the  random  numbers. 

4-  Generates  the  time  axis. 

1-  Calculates  the  FFT  of  the  input  test  signal. 

2-  Calculates  the  FFT  Theoretical  for  comparison. 

1-  Computing  the  filtering  process  in  the 
frequency  domain. 


COD,  B.m 


1-  Calculates  the  IFFT  of  the  filtered  signal. 
1-  Coherent  detection  of  the  filtered  signal. 


C-  THE  PROGRAM 


In  this  section  a  detailed  explanation  of  each  block  is  presented 
for  a  better  understanding  of  how  the  computer  program  works. 


Figure  7  -  Program  structure 


2  1 


1  -  BLOCK  INP_B.m 


The  input  data  block  is  for  introducing  the  variables  needed  to 
run  the  program  like  the  carrier  frequency  (fc),  the  amplitude  of  the 
signal  (Ac),  the  central  frequency  of  the  interference  signal  (fi),  the 
step  interval  in  order  to  simulate  the  movement  of  the  detector 
array  (delta),  the  number  of  loops  that  is  needed  for  the  program  (I), 
and  the  seed  number  for  the  random  number  generator  (seed).  All 
values  are  presented  on  the  screen. 

2-  BLOCK  M_BPSK.m 

This  block  controls  the  sequence  during  the  simulation  and  the 
loop  to  simulate  the  movement  of  the  excisor-detector  array.  In 
order  to  reduce  the  amount  of  memory,  the  data  is  stored  in 
temporarily  files  to  be  used  by  the  other  section.  After  each  block 
has  been  run  and  the  data  has  been  stored  in  a  temporary  file,  the 
memory  is  cleared  before  executing  the  next  block  to  keep  the 
computer  from  running  out  of  memory. 

3  -  BLOCK  FILTER.m 

The  main  functions  of  the  Filter. m  block  are  to  build  the  filter 
transfer  function  and  to  create  the  frequency  axis  that  will  be  used 
by  the  simulation.  The  filter  consists  of  a  signal  that  represents  a 
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detector  array  in  the  frequency  domain.  An  ideal  representation  of 
the  excisor-detector  array  is  showed  in  Figure  8.  In  Figure  9,  the 
filter  magnitude  response  of  the  model  of  the  excisor-detector  array 
is  presented.  The  width  (w)  may  vary  from  0  to  120  KHz  and  the 
depth  (d)  from  0  to  -90  dB.  Figure  8  also  shows  the  center  frequency 
for  each  notch  in  the  excisor-detector  array  assuming  that  the  array 
is  centered  at  70  MHz.  Appendix  C  provides  useful  information  to 
work  with  the  filter,  like  the  frequency  at  which  any  notch  starts  and 
ends  for  a  given  width  and  the  correct  value  of  notch  depth  that 
represents  the  desired  attenuation. 

Performing  this  RF  filtering  process  in  the  frequency  plane  of  an 
optical  spectrum  analyzer  has  advantages  over  the  electronic 
techniques  [Ref.  6].  When  an  object  is  placed  in  the  frequency  plane 
of  the  optical  spectrum  analyzer,  the  magnitude  of  the  affected 
frequency  interval  is  attenuated  and  the  phase  remains  constant. 
Also,  many  filters  may  be  realized  simply  by  varying  the  size  of  the 
notch  in  the  frequency  plane. 

The  excisor-detector  is  placed  at  one  focal  length  away  from  the 
Fourier  transform  lens  into  the  frequency  plane.  It  has  13 
photodiodes  elements  that  represent  the  different  channels  750  KHz 
wide,  adding  to  a  total  9.75  MHz  window  as  shown  in  Figure  8. 
(These  values  represent  values  used  in  an  actual  laboratory  setup). 
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Figure  9  -  Filler  magnitude  response 

In  the  computer  program  (Appendix  B,  p.  66)  the  notch  depth  is 
selected  to  represent  the  desired  attenuation.  The  separation 
between  points  in  the  frequency  domain  is  5  KHz.  To  represent  a 
notch  60  KHz  wide,  12  numbers  or  different  elements  are  needed. 

The  frequency  axis  is  created  considering  that  the  sampling 
frequency  (fs)  must  be  greater  than  twice  the  highest  frequency 
present  in  the  signal  that  is  being  sampled.  That  is,  the  sampling 
frequency  must  satisfy  the  inequality  fs>2fmax  [Ref.  7].  The  highest 
frequency  in  our  test  signal  is  fmax=70  MHz,  so  the  sampling 
frequency  must  be  greater  than  the  Nyquist  frequency  2fmax=140 
MHz.  For  the  simulation  fs=  163.84  MHz  has  been  chosen. 


Using  MATLAB386  [Ref.  4]  the  following  expression  was  used  to 
create  the  frequency  axis: 

f  =  (bin_number  -l)*fs/N  (10) 

where 

f  =  frequency  axis, 
fs  =  sampling  frequency, 

N  =  number  of  points  or  vector  length,  and 
bin_number  =  N/2+1. 

In  order  to  represent  the  filter  signal  centered  at  70  MHz  the 
values  used  in  the  previous  expression  have  been  N  =  32768, 
bin_number-l  =  16384,  and  fs  =  163.84  MHz. 

Two  plots  are  presented  on  the  screen,  showing  the  detector 
array  filler  signal  and  the  filter  magnitude  response. 

4-  BLOCK  S.BPSK.m 

This  block  constructs  the  input  test  BPSK  signal,  the  interference 
signal,  and  the  time  axis  to  be  used  in  the  program. 

The  BPSK  signal  is  represented  by  [Ref.  7] 

s(t)  =  Ac  cos[(Oc  t+Dp  m(t)]  (11) 

where  m(t)  is  a  bipolar  baseband  data  signal  with  peak  values  ±1. 
Recalling  that  cos(x)  and  sin(x)  are  even  and  odd  functions  of  x,  the 
BPSK  signal  can  be  obtain  by  [Ref.  7] 
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s(t)  =  (Ac  cos  Dp)  cos  COc  t  -  (Ac  sin  Dp  )  m(t)  sin  (Oc  t  (12) 

first  term  second  term 

where  the  first  term  is  called  the  pilot  carrier  term  and  its  level  is 
set  by  the  value  of  the  peak  deviation  A0  =  Dp, 

For  digital  angle-modulated  signals,  the  digital  modulation  index, 
h,  is  defined  by  [Ref.  7] 

h=2Ae/7r  (13) 

where  2  A0  is  the  maximum  peak-to-peak  phase  deviation  during 
the  time  that  it  takes  to  send  one  symbol,  Ts.  For  binary  signaling, 
the  symbol  tim-"  equal  to  the  bit  time,  Ts  =  Tb. 

The  po*v  ,,  in  the  second  term  (called  the  data  term)  needs  to  be 
maximized.  This  is  accomplished  by  letting  A0  =  Dp  =  90°  [Ref.  7], 
which  corresponds  to  a  digital  modulation  index  of  h  =  1,  For  this 
case,  changing  sin  by  cos,  the  BPSK  signal  becomes 

s(t)  =  Ac  m(t)  cos  (Oct  (14) 

where 

(Oc  =  2  7C  fc,  and 
Ac  =  Amplitude. 

Equation  14  is  the  algorithm  applied  in  the  test  signal 
generation.  The  m(t)  signal  is  a  random  binary  data  that  represents 
the  modulating  waveform  originated  from  a  digital  information 
source  and  is  illustrated  in  Figure  10.  The  use  of  random  data  should 
be  more  effective  than  using  any  other  deterministic  method. 
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Figure  10  -  Binary  random  data 

The  random  number  generator  (RNG)  employed  must  fulfill  the 
following  requirements  [Ref,  8]: 

•  Uniform:  It  must  pass  several  statistical  test  for  U(0,1)  uniform 

distribution. 

•  Reproducible:  To  be  useful,  the  RNG  should  reproduce  the 

same  string  of  numbers  given  the  same  seed  every  time;  this 
is  of  particular  importance  for  testing  purpose. 

•  Minimum  Memory:  It  should  be  efficient  to  use  and  fast. 

•  Non-repeating:  A  number  sequence  must  not  repeat  itself 

during  a  simulation  run. 

•  Statistical  Independent;  This  property  is  very  difficult  to 

achieve  because  of  the  recursive  properties  of  the  generator. 
The  mathematical  algorithm  was  based  on  the  linear 

congruential  method  of  the  form  (Ref.  8]: 
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zi  =  (a*zi-i  +  c)  MOD  m 


(15) 


where 

a  =  multiplier, 
c  =  increment, 
m  =  modulus,  and 
zo  =  seed. 

It  can  be  shown  that  the  pseudorandom  number  generator  can 
be  constructed  using  the  following  expression  [Ref.  8] 

zi=  (16807  *zi-i)  MOD  (16) 

For  the  random  case  the  test  signal  generated  by  Eq.  14  will 
produce  a  continuous  spectrum  as  will  be  presented  later. 

The  interference  signal  is  also  generated  in  the  present  block 
based  on  a  sinusoidal  signal  with  an  amplitude  three  times  the 
original  BPSK  signal  (in  the  time  domain)  and  a  carrier  frequency 
selectable  from  the  input  data.  As  a  result  of  adding  those  two 
signals,  the  test  BPSK  signal  plus  interference  is  produced. 

The  generation  of  the  time  axis  is  other  function  of  this  block.  A 
total  time  of  0.2  milliseconds  with  intervals  of  6.1  nanoseconds  has 
been  selected. 

Plots  presenting  the  random  binary  data,  BPSK  test  signal,  and 
BPSK  test  signal  plus  interference  are  presented  on  the  screen  as 
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shown  in  Figures  11  and  12.  A  carrier  frequency  of  70  MHz  and  an 
amplitude  of  1  volt  was  used  for  the  BPSK  signal  and  a  central 
frequency  of  67.4  MHz  and  3  volts  amplitude  for  the  interference 
signal. 


Figure  1 1  -  BPSK  signal  and  modulated  waveform 

5  -  BLOCK  FFT.B.m 

The  FFT_b.m  block  uses  the  input  signal  constructed  in  the 
previous  block  and  calculates  the  fast  Fourier  transform.  The  first 
step  in  applying  the  discrete  Fourier  transform  is  to  choose  the 
number  of  samples  N  and  the  sampling  interval  T,  and,  secondly,  to 
interpret  the  results  correctly.  The  discrete  transform  approximates 
the  continuous  transform  rather  poorly  for  the  higher  frequencies 
[Ref.  7].  To  reduce  this  error  it  is  necessary  to  decrease  the  sample 
interval  T  and/or  to  increase  N. 
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(17) 


In  our  case  the  relationship  used  was 

T  =  Period  /  N 
where 


N  =  32768  (  N  =  2*^  where  (p  is  an  integer  value),  and 
Period  =  0.2  millisec. 

So  T  =  6.1  nanosec. 
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Figure  12  -  BPSK  signal  plus  interference 


From  MATLAB386  [Ref.  4]  the  FFT(x)  is  the  discrete  Fourier 
transform  (DFT)  of  vector  x,  computed  with  a  radix-2  fast  Fourier 
transform  algorithm.  If  the  length  of  x  is  not  an  exact  power  of  two 
then  it  is  padded  with  trailing  zeros.  This  is  not  our  case  because  we 
have  selected  N  =  32768  =  2^^. 

The  function  that  implements  the  transform  is  given  by 

N-  1 

XCk-^l)=  X  x(n  -H) 

n  =0 


3  1 


(18) 


where  Wj^  =  and  N  =  length(X).  Note  that  the  series  (18)  is 

written  in  an  unorthodox  way,  with  the  subscript  from  k+1  instead  of 
the  usual  k  because  subscripts  on  Matlab  vectors  run  from  1  to  N 
instead  of  from  0  to  N-1  [Ref.  4]. 

In  order  to  plot  the  FFT  results,  the  frequency  axis  generated  in 
the  Filter. m  block  is  used.  Actually  the  power  spectral  density  (PSD) 
is  calculated  and  plotted  as  shown  in  Figures  13. 

Also,  for  comparison  purposes,  the  theoretical  power  spectral 
density  of  a  BPSK  signal  with  the  same  characteristics  as  the  test 
signal  was  calculated  and  plotted  using  the  following  algorithm: 


P(f)  = 


A  fsinWf-gT^) 


4R|  ll(f-gT^ 


where 


Tb  =  bit  time  =  1/R,  and 
R  =  bit  rate. 


The  resulting  BPSK  spectrum  is  shown  in  Figure  14.  The  null-to- 
null  bandwidth  is  2R.  The  figure  shows  plots  of  the  calculated  and 
the  theoretical  PSD  of  the  test  signal.  As  we  can  see,  the  calculated 
spectrum  (solid  line)  appears  to  be  noisy  but  follow''  the  theoretical 
spectrum  plotted  using  equation  19  (dotted  line).  Therefore  the 
algorithm  used  is  a  good  representative  of  a  BPSK  signal. 
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Figure  13  -  PSD  of  a  test  BPSK  signal  plus  interference 


Figure  14  -  Theoretical  and  calculated  PSD  of  a  BPSK  signal 


6  -  BLOCK  S_FlL.m 


The  purpose  of  this  block  is  to  take  the  FFT  of  the  input  signal 
and  multiply  it  with  the  detector  array  filter.  As  stated 


before,  the  photodiode  array  with  the  small  space  between  each  of 
the  photodiodes  can  produce  the  filter.  Placing  a  notch  where  the 
interference  is  present  produces  a  narrowband  frequency  rejection 
via  a  notch  filter.  Also  the  harmful  effect  of  notches  on  the  other 
frequency  components  of  the  original  signal  must  be  considered. 

A  plot  of  the  filter  transfer  function  and  the  test  signal  spectrum 
will  be  presented  in  order  to  visualize  the  effect  of  the  filter  on  the 
signal. 

Finally  a  comparison  between  the  filtered  signal  spectrum  and 
the  original  signal  spectrum  is  showed  to  evaluate  a  priori  the 
filtering  effect. 


rpEQUENCy,  Hi  ilO'’ 

Figure  15  -  Filler  transfer  function  and  lest  signal  spectrum 

An  example  of  these  plots  are  illustrated  in  Figures  15  and  16. 
In  Figure  15  the  interference  dominates  the  spectrum  and  the  BPSK 
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signal  appears  relativelly  small.  Also  the  magnitude  of  the  filter 
response  of  the  excisor-detector  array  is  shown.  The  overlap  of  the 
interference  with  the  third  notch  is  clearly  shown.  After  the  filtering 
process  takes  place  (Figure  16b),  the  BPSK  signal  increases  with 
respect  to  the  interference  signal,  showing  the  effectiveness  of  the 
filtering  process. 


b  -  Filtered  signal  spectrum. 

Figure  16  -  Filtered  and  original  signal  spectrum  comparison 


7  -  BLOCK  IFFT_B.m 

Here  the  output  signal  is  determined  by  an  inverse  transform 
of  the  filtered  signal  calculated  in  the  previous  block.  In  order  to 
perform  the  inverse  fast  Fourier  transform  (IFFT)  the  following 
algorithm  was  used  [Ref.  4]; 


1  '  -k 

x(n  +1)=-^  Z  X(k+1) 


The  recovered  signal  is  presented  and  compared  with  the 
original  test  BPSK  signal  without  interference  as  shown  in  Figure  17, 


ORIGINAL  DPSK 


0  OZ  Oi  06  00  I  I.Z  1.4  16  10  Z 

TIMR.  sec  slO"^ 

Figure  17  -  Original  and  recovered  BPSK  signal 

We  can  see  that  the  recovered  signal  experiences  a  distortion 
due  to  the  notches  and  the  windowing  effect  of  the  excisor-detector 
array,  in  addition  to  the  interference  signal. 
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8  -  BLOCK  COD_B.m 


The  optimum  detection  of  a  BPSK  signal  consists  of  coherent 
detection  as  showed  in  Figure  18.  The  gated  integrator  that  is  turned 
on  at  the  beginning  of  each  bit  period  (Tb)  and  integrates  to  the  end 
of  the  bit.  Then  the  integrator  output  is  sampled  and  mapped  into 
digital  values  (a  binary  1  if  the  sample  value  is  a  positive  voltage 
and  a  binary  0  if  the  sample  value  is  negative).  The  integrator  is 
then  reset  to  a  zero  initial  condition  and  this  procedure  is  repeated 
for  detection  of  the  next  bit.  This  optimum  detector  is  called  an 
integrate  and  dump  detector  [Ref.  7)  presented  in  Figure  19. 


Figure  18  -  Coherent  detection 


The  operation  described  is  simulated  with  a  FOR  loop  that 
computes  the  cumulative  sum  during  the  bit  time  and  is  dumped  to 
zero  at  the  end.  To  map  the  integrator  output  (Figure  20)  into 
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digitals  values,  another  loop  with  logic  statements  is  performed, 
reading  the  middle  point  of  each  bit  time.  These  operations  generate 
a  string  of  binary  data  that  will  be  compared  to  the  original  random 
binary  information  as  shown  in  Figure  21. 


Figure  19  -  Integrate  and  dump  detector 


TIME,  sec 


iio-< 


Figure  20  -  Inlcgralor  oulpul  (Tb  =  0.78  jiscc) 
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0  0  2  O'*  0  G  0  0  1  1.2  14  16  in  2 

TIME,  9tc  Il0*< 

Figure  21  -  Original  and  recovered  binary  data. 

In  order  to  detect  if  an  error  is  present,  the  original  random 
binary  data  and  the  recovered  binary  data  are  compared.  This 

information  is  presented  on  the  screen  showing  a  1  if  no  error  is 
detected  otherwise  a  -I  represents  an  error  (Figure  22). 


NO  EfiROR 


%  -0.5 
< 


f>  0  2  0  1  0  r>  on  i  12  \a  i  o  in  2 


Figure  22  -  Error  detection 
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If  in  the  input  data  the  number  of  intervals  (I)  selected  is 
greater  than  1,  the  program  regenerates  a  test  BPSK  signal  with  the 
interference  frequency  shifted  an  amount  equal  to  the  frequency 
interval  (delta)  introduced  at  the  lNP_B.m  block.  This  action 
simulates  the  movement  of  the  excisor-detector  array  in  order  to 
evaluate  the  effect  of  this  translation  on  the  output  signal. 
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V.  PROGRAM  EVALUATION 


An  evaluation  of  the  computer  program  presented  in  Chapter  IV 
and  Appendix  B  was  performed  under  different  conditions  as  shown 
in  Table  II.  The  main  idea  was  to  be  sure  that  the  output  for  each 
test  agreed  with  what  we  expected  or  with  previous  experimental 
results. 


TABLE  II  -  TEST  CONDITIONS 


Excisor- 

^\detector 

Test 

Signal 

YES 

with  notches 

NO 

without  notches 

BPSK  signal 
without 
interference 

Degradation  of 
the  original 
signal  due  to 
notches 

Degradation  of 
the  signal 
passing  through 
the  window 

BPSK  signal 
with 

interference 

Filtering 
capability  over 
the  interference 

Degradation  of 
the  signal. 

Errors  on  the 

Data  recovered 

A-  TEST  N“  1 


For  the  purpose  of  the  test,  a  BPSK  signal  without  interference 
was  consider.  An  excisor-detector  array  without  notches  (ideal,  as  a 
window  9.75  MHz  wide  using  an  attenuation  factor  of  1  (0  dB))  was 


used  to  represent  the  detector  array  centered  at  70  MHz.  The  BPSK 
signal  generated  as  a  test  signal  had  the  following  parameters  that 
remained  constant  during  all  other  tests. 


TABLE  III  -  TEST  SIGNAL  PARAMETERS 


Parameter 

Description 

Value 

A 

Amplitude  BPSK  signal 

1  volt 

Ai 

Amplitude  interference 

3  volts 

Tb 

Bit  Time 

0.39  nsec 

fc 

Carrier  Freq.  BPSK 

70.0  MHz 

fi 

Center  Freq.  Interference 

67.4  MHz 

seed 

PRNG  seed  number 

236578 

In  the  first  test  the  central  frequency  of  the  interference  was 
selected  as  fi=0  in  order  to  simulate  a  no-interference.  The  value  of 
the  bit  time  (Tb)  equal  to  0.39  psec  was  employed  to  put  the  minimal 
amount  of  the  original  signal  on  the  window  (main  lobe  and  the  first 
side  lobes).  The  results  of  this  simulation  are  summarized  below 
(Figure  23)  and  match  what  was  expected. 


AMPUTUDE.  volt,  VMPLlTt  DE.  volt,  pO^j-R 


a-  PSD  of  ihc  Test 
signal.  The 
ihcorclical  PSD  is 
shown  in  dolled 
line  and  the 
calculated  in  solid 
line. 
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b-  Filler  (cxcisor- 
deteclor  array) 
signal  and  the 
spectrum  of  the 
test  signal. 
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c-  Spectrum  of 
the  test  BPSK 
signal  after 
filtering 
(frequency 
domain). 
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d-  Recovering  the 
BPSr.  .  and 
comparing  it  with 
the  original  test 
BPSK  signal. 
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e-  Integrator 
output  during  the 
detection  process. 
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f-  Plot  of  the 
original  Random 
Binary  Data  m(t) 
and  the  recovered 
data  after  the 
process. 
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g-  Checking  for 
bit  error. 
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Figure  23  -  Output  Test  N-  1 

As  was  expected  the  degradation  of  the  test  BPSK  signal  due  to 
the  window  9.75  MHz  wide  was  very  small  and  the  data  information 
was  recovered  without  any  error. 

B-  TEST  N«  2 

The  second  test  used  the  same  test  BPSK  signal  as  before  with 
the  addition  of  an  interference.  The  excisor-detector  was  still  a 
window  9.75  MHz  wide  with  an  attenuation  factor  of  1  (0  dB).  Figure 
24  shows  the  results  of  Test  N°  2. 
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a-  PSD  of  ihc  icsl 
BPSK  signal  plus 
interference  is 
presented  in  solid 
line.  Theoretical 
plot  is  shown  in 
dotted  line  as  a 
reference.  The 
interference  is  30 
db  above  the  test 
signal . 
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b-  Filter  signal 
and  spectrum  of 
the  test  BPSK 
signal  plus 
interference  is 
presented.  Note 
that  the 
interference 
signal  dominates 
the  spectrum 
when  a  linear 
scale  is  used. 
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c-  Output  after 
passing  the  test 
signal  through 
the  window. 
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e-  Output  of 
integrating 
detector. 
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f-  Comparison  of 
the  original 
random  binary 
data  m(t)  and  the 
recovered  data. 
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g-  Error  detection 
output. 
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Figure  24  -  Output  Test  N®  2 

As  seen  in  Figures  24f  and  24g,  the  interference  signal  produces 
a  significant  distortion  on  the  recovered  BPSK  signal  as  well  as  on  the 
random  binary  data.  Hence  a  filter  is  desired  to  reduce  the  harmful 
interference. 


48 


C-  TEST  N“  3 


From  the  previous  results  an  effective  filtering  process  is 
needed.  The  actual  excisor-detector  array  will  have  notches  in  the 
frequency  domain.  We  want  to  examine  the  effects  of  these  notches 
on  the  BPSK  signal  alone.  The  present  test  was  performed  with  a  test 
BPSK  signal  without  interference  passing  through  an  excisor-detector 
with  notches  60  KHz  wide  and  with  a  depth  of  20  dB  (see  Figure 


25b). 


a-Thc  PSD  of  the 
test  signal 
without 
interference 
within  the 
window  of  the 
excisor-detector 
array. 
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b-  Filter  signal 
and  the  spectrum 
of  the  test  BPSK 
signal  without 
interference. 

Note  the  notches 
over  the 

broadband  signal. 
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c-  After  passing 
through  the 
excisor-dctcctor 
filter  signal  the 
spectrum  of  the 
filtered  signal  is 
presented. 
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d-  Comparison 
between  the 
original  and 
recovered  BPSK 
test  signal. 
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e-  Integrator 
output  during  the 
detection  process. 


No  significant  distortion  was  observed  between  the  results  of 
this  filtered  signal  and  the  unfiltered  BPSK  signal  (Test  N°  I).  It  is 
reasonable  to  state  that  the  effect  of  these  notches  (60  KHz  wide  and 
20  dB  depth)  is  negligible  on  the  broadband  signal. 


D-  TEST  N»  4 


Finally  we  consider  a  common  situation.  The  test  BPSK  signal 
was  considered  with  an  interference  present  and  the  excisor-detector 
array  had  notches  between  detectors  that  will  produce  the  filtering 
effect.  A  notch  was  aligned  on  the  interference  frequency  (see  Figure 
26b). 


a-  PSD  of  a  test 
BPSK  signal  with 
an  interference 
at  67.4  MHz  (solid 
line)  and  the 
theoretical  PSD  of 
a  BPSK  signal  as  a 
reference . 
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c-  Spectrum  of 
the  filtered 
signal.  The 
interference  is 
attenuated  and 
the  signal 
spectrum 
increases, 
relative  to  the 
interference. 


d-  The  recovered 
filtered  BPSK  and 
the  original 
signal  without 
interference. 


c-  The  integrator 
detector  output. 
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f-  The  original 
data  and  the 
recovered  data. 
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g-  Error  detection 
between 
recovered  data 
and  test  data  is 
performed . 
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Figure  26  -  Output  Test  N®  4 


During  this  test  the  effect  of  the  narrowband  notches  was 
evaluated  showing  that  this  filtering  technique  was  effective  against 
the  narrowband  interference  signal  considered. 
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E-  OTHER  TESTS 


1  -  NO-FILTER  TEST 

An  early  test  was  performed  to  check  the  computer  program. 
The  test  signal  was  the  BPSK  signal  without  interference.  No  filtering 
process  was  included.  We  only  generated  and  recovered  the  BPSK 
signal  (Appendix  B,  p.  66).  The  signal  was  recovered  without  any 
distortion  (Figs  27,  28  and  29). 
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Figure  27  -  Original  and  recovered  BPSK  signal 
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Figure  28  -  Original  and  recovered  binary  data 
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Figure  29  -  Error  detection  output 


2  -  MOVING  THE  EXCISOR-DETECTOR 


One  of  the  objectives  was  to  simulate  motion  of  the  detector 
array  to  study  the  effects  of  the  notch  passing  over  the  interference 
signal.  For  this  situation  the  program  was  provided  with  a  loop  that 


simulated  this  movement.  On  the  following  figures  the  filtering 
effect  can  be  seen.  As  the  interference  passes  through  the  notch  (60 
KHz  wide  and  20  dB  depth),  different  attenuation  effects  take  place. 
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Figure  30  -  Test  signal  spectrum  (notch  at  67,35  MHz) 
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Figure  31  -  Test  signal  spectrum  (notch  at  67.38  MHz) 
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Figure  34  -  Test  signal  spectrum  (notch  at  67.41  MHz) 
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Figure  35  -  Test  signal  spectrum  (notch  at  67.42  MH^) 
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As  we  can  see  in  Figures  30  to  35,  a  good  filtering  effect  is 
obtained  using  a  relatively  narrow  notch  with  an  reasonable 
attenuation.  The  interference  signal  that  initially  dominates  the 
spectrum  was  progressively  attenuated  as  the  notch  was  moved  over 
it  (see  Figure  32). 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  has  presented  the  basic  principles,  characteristics  and 
potential  capabilities  of  optical  excisors.  The  program  developed 
using  MATLAB386  provides  a  consistent  tool  to  aid  in  the  study  of 
acousto-optic  spectrum  analyzer  and  excision  techniques. 

Processing  broadband  signals  requires  the  removal  of  strong, 
narrowband  interfering  signals.  Narrowband  interference  has  been 
eliminated  using  a  Bragg  cell  receiver  with  an  excisor-detector  at  the 
frequency  plane.  Moving  the  excisor-detector  by  optical  and 
mechanical  means  made  it  possible  to  place  a  notch  with  a  given 
characteristic  over  the  narrowband  interference  signal.  The 
temporal  characteristics  of  the  broadband  signal  are  recovered  by 
imaging  the  filtered  light  distribution  and  a  reference  beam  onto  a 
photodetector  array  to  allow  coherent  detection. 

All  the  results  produced  with  the  program  during  different  tests 
agreed  with  the  values  obtained  from  available  references. 

Areas  for  potential  use  of  the  program  developed  are  the 
following: 

1.  Evaluation  of  the  excision  effect  varying  the  notch  width  and 
depth  using  the  bit  error  rate  as  an  indicator  of  its  effectiveness. 

2.  Testing  different  input  signals  and  narrowband  interference 
signals. 
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APPENDIX  A  -  NOTATION  USED  IN  THE  PROGRAM 


The  notation  used  consistently  throughout  this  thesis  is 
summarized  for  the  mayor  abbreviations  and  symbols. 


A  = 

Ac  = 
b  = 

B  = 
c  = 

delta  = 
DET  = 
f  = 
fc  = 
ft  = 

F  = 

Fdb  = 
HL  = 

G  = 

G1  = 

H  = 

I  = 
k  = 
m  = 
mo  = 
n  = 

0  = 

Psfl  = 
Pssl  = 
PSS  = 
PS2  = 


Total  time 

Binary  random  voltage  (Vector) 

Signal  amplitude 
Number  of  binary  data 
Integrator  Output  (Vector) 

Increment  (PRNG) 

Frequency  interval 
Coherent  detection  (Vector) 

Frequency  axis  (Vector) 

Carrier  frequency 
Interference  central  frequency 
Channel  (Vector) 

Detector  array  filter  signal  in  db  (Vector) 

Detector  array  filter  signal  (Vector) 

PSD  of  a  BPSK  signal  theoretical  norm. (Vector) 

PSD  of  BPSK  signal  theoretical  in  dB  (Vector) 
Recovered  signal  mapped  into  binary  data  (Vector) 
Number  of  loops  in  the  program 
Generic  counter 

Number  of  samples  per  bit  time  (Tb) 

Modulus  (PRNG) 

Total  number  of  elements  in  a  vector 
Right  half  of  a  notch  (Vector) 

PSD  of  filtered  signal  normalized  (Vector) 

PSD  of  signal  plus  interference  normalized  (Vector) 
PSD  of  signal  plus  interference  in  dB  (Vector) 

PSD  of  signal  without  interference  normalized  (Vec.) 


PSSl  = 

Q  = 

rs  = 

RN  = 
s  = 
sa  = 
si  = 
sd  = 
seed  = 
S  = 

SD  = 
SF  = 

t  = 
u  = 
zi  = 


PSD  of  signal  without  interference  in  dB  (Vector) 
Left  half  of  a  notch  (Vector) 

Recovered  signal  in  time  domain  (Vector) 
Random  number  (0  -  1) 

BPSK  signal  (Vector) 

Integrator  and  dump 
Interference  signal  (Vector) 

Signal  (BPSK)  plus  interference  (Vector) 

Seed  for  the  pseudorandom  number  generator 
FFT  of  BPSK  signal  without  interference  (Vector) 
FFT  of  BPSK  signal  plus  Interference  (Vector) 
Signal  filtered  (Vector) 

Time  axis  (Vector) 

Generic  counter 

Remainder  after  division  (PRNG) 
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APPENDIX  B  -  COMPUTER  PROGRAM 


The  following  are  the  different  blocks  that  using  MATLAB386  as 
a  master  program  make  possible  to  perform  the  simulation  and 
evaluate  the  results.  The  titles  of  each  block  in  bold  typeface  are  not 
part  of  the  program. 


BLOCK  Input  initial  Data  (INPJB.m) 


echo  on 
clc 

%  THIS  IS  A  PROGRAM  DEVELOPED  IN  ORDER  TO  EVALUATE 
%  THE  FILTERING  CAPABILITY  OF  A  NOTCH  PRODUCED  BETWEEN 
%  CHANNELS  IN  A  DETECTOR  ARRAY. 

%  NOW  YOU  WILL  BE  REQUESTED  TO  INTRODUCE  SOME  PARAMETERS 


pause  9c  Strike  any  key  to  continue 

clc 

echo  on 

%  THE  CARRIER  FREQUENCY  MUST  BE  SELECTED  BETWEEN  65  AND  75  MHz 
%  WE  SUGGEST  70  MHz  IN  ORDER  TO  WORK  IN  THE  CENTER  OF  THE 
%  DETECTOR  ARRAY. 

echo  off 

inputf'Which  is  the  Carrier  Frequency  (Hz)?  '); 
fc=ans;  %  Carrier  frequency 

clc 

echo  on 

%  THE  AMPLITUDE  OF  THE  SIGNAL  MAY  BE  VARIED 
%  WE  SUGGEST  TO  USE  1  volt. 

echo  off 

inputCSelect  the  Amplitude  of  the  Signal  =  ’); 

Ac=ans;  %  Signal  amplitude 


clc 

echo  on 

%  THE  FREQUENCY  OF  THE  INTERFERENCE  SIGNAL  MUST  BE  IN  THE 
%  FREQUENCY  RANGE  OF  THE  DETECTOR  ARRAY  OTHERWISE  WILL  BE 
%  AUTOMATICALLY  FILTERED  (65  TO  75  MHz). 
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echo  off 

inputCSelect  the  Interference  Frequency  (Hz)  =  '); 

fi=ans;  %  Interference  Frequency 

clc 

echo  on 

%  IN  ORDER  TO  RUN  THE  PROGRAM  AND  TO  SIMULATE  THE  MOVEMENT 
%  OF  THE  DETECTOR  ARRAY  YOU  HAVE  TO  SELECT  THE  STEP  IN 
%  FREQUENCY  THAT  YOU  WANT  TO  CONSIDER  FOR  EACH  LOOP 
%  FOR  EXAMPLE:  YOU  WANT  TO  MOVE  THE  DETECTOR  ARRAY  IN 
%  STEPS  OF  10  KHz.  INSERT  THIS  VALUE  IN  Hz  (.10000  ) 

echo  off 

inputCSelect  the  frequency  interval  (Hz)  =  '); 
delta=ans;  %  Frequency  Interval 

input('How  many  intervals  do  you  want?  '); 

I=ans;  %  Number  of  loops 


clc 

echo  on 

%  TO  PROVIDE  THE  PROGRAM  WITH  ONE  ESSENTIAL  CHARACTERISTIC 
%  IN  TEST  AND  EVALUATION  (IT  MUST  BE  REPRODUCIBLE)  IS  NECESSARY  % 
TO  INTRODUCE  A  SEED  FOR  THE  PSEUDORANDOM  NUMBER  GENERATOR 
(PRNG). 

echo  off 

inputCInput  a  SEED  for  the  Random  Number  Generator  =  '); 

seed=ans;  %  SEED  for  PRNG 

echo  on 

save  tcmp_INPUT  fc  fi  Ac  delta  I 

M_BPSK  %  Calling  the  main  control  block 

end 


Bt  OCK  Main  Control  (M_BPSK.m) 


echo  off 

clear 

FILTER 

clear 

clc 

echo  off 

load  temp_INPUT 

for  k=0;I 

load  temp_INPUT 
fc=fc+delta*k 

fi  =  fi  +  delia*k 
pausc(3 ) 


%  Turns  program  display  off 
%  Clears  memory 
%  Develop  filter  signal 


%  Loop  for  frequency  steps 

%  Carrier  frequency  plus 
%  interval 

%  Interference  frequency 
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S_BPSK 

% 

Generate  the  input  signal 

clear 

% 

Clears  memory 

FFT_B 

% 

Develops  FFT  of  input 

% 

signal 

clear 

% 

Clears  memory 

S_FIL 

% 

Calculate  the  resulting 

% 

signal  passing  through  the 

% 

detector  array 

clear 

% 

Clears  memory 

IFFT_B 

% 

Develops  the  IFFT  of  the 

% 

filtered  signal 

clear 

% 

Clears  memory 

<X»_B 

% 

Coherent  detection 

clear 

% 

Clears  memory 

end 

% 

End  loop 

end 

% 

End  of  main  program 

cla 

% 

Clears  program  from  memory 

BLOCK  Detector  Array  Filter  Signal  (FILTER. m) 


load  temp_S 

axis; 

cic 


%  Is  time  to  create  the  detector  array  filter  signal 

%FIL=[ones(l ,32768)];  %  Simulates  no-filter 

Q=[l  1  1  1  1  1  1  1  1  1  1  Ij; 

%  left  1/2  notch 
0=[1  1  1  1  1  1  1  1  1  1  1  Ij; 

%  right  1/2  notch 

F=[Q  ones(l,126)  Oj; 

FlL=[zeros(l,  13028)  0:1:1  ones(1.137)  .. 
OFFFFFFFFFFF.. 

Q  ones(l,l37)  1:-1:0  zeros/ 1,..  %  Detector  array 

2808)  0:1:1  ones(l,137)  ..  %  Filler  signal 

OFFFFFFFFFFF.. 

Q  ones/ 1,1 37)  1:-1:0  zeros/ 1,1 3028)]; 
pause/3) 

c  Ic 

%  We  need  to  create  a  frequency  axis 
f=163.84e6*/0: 16383)/32768;  %  Frequency  axis 

clc 

%  The  following  is  a  plot  of  the  detector  array 
%  Filter  Signal 

axis/[65c6  75e6  0  1.2!); 
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pIot(f.FIL(I:I6384)),.. 

tillcCDETECTOR  ARRAY  RLTER  SIGNAL’)... 
xlabeK'FREQUENCY.  Hz')... 
ylabeK’AMPLlTUDE.  volts'). .. 
pause(5) 

clg 

clc 

%  Now  we  look  at  the  filter  magnitude  response 

FILl=FIL+eps: 

FdB=10*logl0(FILl); 

axis((64€6  76e6  -70  10]); 

ploi(f.FdB(l:16384))... 

titleCFILTER  MAGNITUDE  RESPONSE’)... 

xlabeK'FREQUENCY.  Hz'),.. 

ylabel('lH(f)l,  dB')... 

pause(5) 

clg 

clc 

save  temp_F  FIL  f 
end 


BLOCK  BPSK  Signal  Generator  (S_BPSK.m) 


echo  off 
load  temp_F 
clc 

a=0.0002;  %  Total  ti.ne 

t=0:a/32767.5;a;  %  Time  axis 

echo  on 
clc 

%  In  order  to  produce  a  BPSK  signal  we  need  to  generate 
%  a  binary  random  voltage  with  values  +1  or  -1. 


c=16807; 

% 

Multiplier 

mo=2''31-l; 

% 

Modulus 

zi=seed; 

% 

Remainder  after  division 

n=32768; 

% 

Total  #  of  elements  in  the  vector 

m=64; 

% 

#  of  samples  per  bit  time 

b=512; 

% 

#  of  binary  data 

A=ones(l  ,n); 

% 

Base  vector 

p-u, 

for  k=m;m:n 

% 

Loop 

g=rem(c*zi,mo); 

% 

Pseudorandom  number 

RN=g/mo; 

% 

generator 
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zi— gi 

d=(2*round(RN))-l;  %  [  +1  or  -1] 

A(l.(p*ni)+l:k)=A(l.(p*m)+l;k)*d; 
p=p+l: 

end  %  End  of  Loop 

%  Now  is  lime  to  plot  the  binary  data 

axis([0  a  -Ac-1  Ac+1]); 
plot(t,A),.. 

tilleCBINARY  RANDOM  DATA'),.. 
xlabeK'TIME.  sec')... 
ylabelCAMPLITUDE.  volts')... 
pause 

clc 

%  Next  we  can  form  a  BPSK  Signal  with  a  carrier  frequency  (fc) 

s=Ac'‘‘ A.*cos(2*pi*fc*t);  %  BPSK  Signal 

si=3*Ac*sin(2*pi*l*fi).  %  Interference 

sd=s+si:  BPSK  signal  plus  Interference 

clc 

%  Now  we  are  ready  to  plot  the  pulse  and  the  BPSK  signal 

axis({0  0.00002  -Ac-1  Ac+1)); 
plot(t.s(l:32768).l.A(l:32768).'-.')... 
tiilef'BPSK  TEST  SIGNAL')... 
xlabel('TIME.sec')... 
ylabeU'AMPLlTUDE,  volts')... 
pause 

clg 

clc 

%  Plot  the  BPSK  signal  plus  interference 
axis((0  0.00002  -Ac-1  Ac+1)); 

plot(t.sd).tille('BPSK  SIGNAL  PLUS  INTERFERENCE')... 

xlabelf 'TIME,  sec'),.. 

ylabel('AMPLITUDE')... 

pause 

clg 

clc 

save  temp_S  m  t  s  Ac  fc  a  sd  A  n  b 
save  lemp_SR  m  t  Ac  A  n  b  a 
end 
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BLOCK  FFT  of  BPSK  Signal  (FFT_B.m) 


echo  on 
load  temp_S 
load  temp_F 

clc 

%  Now  we  gel  the  Fourier  transform  of  the  input  signal  (BPSK) 

S=fft(s);  %  FFT  BPSK  Signal 

SD=fft(sd);  %  FFT  BPSK  Signal  +  Interference 

clc 

%  In  order  to  plot  the  power  spectral  density  of  the  signal  we  need  to 
%  make  some  calculations 

Pss=SD.*conj(SD); 

Pssl=Pss/max(Pss);  %  PSD  Signal  +  Interference 

PSS=10*logl0(Pssl);  %  PSD  in  dB 

Psl=S.*conj(S); 

Ps2=Psl/max(Psl );  %  PSD  Signal  without  interference 

PSSl  =  10*logl0(Ps2);  %  PSD  in  dB 

clc 

%  Plot  the  power  spectral  density  of  the  BPSK  signal 
%  plus  interference 

axis([fc-5e6  fc+5e6  -90  0]); 
plot(f,PSS(l:16384)),.. 

lilleCPSD  OF  BPSK  SIGNAL  PLUS  INTERFERENCE’),.. 

xlabcK'FREQUENCY.  Hz')... 

ylabelCPOWER,  dB'),.. 

pause 

hold  on 

clc 

%  For  comparison  we  calculate  the  theoretical  PSD  of  a  BPSK  signal 

R=(l/((a/32768)*m)); 

L=(pi*(f- fc)/R)-f-eps; 

K=(Ac^2)/(4*R); 

X=(sin(L)./L).^2; 

PSD=K*X; 

G=PSD./max(PSD); 

Gl  =  10*logl0(G); 

clc 
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%  Now  we  can  plot  the  theoretical  PSD  of  the  BPSK 


plot(f,Gl(l:16384))... 
meta  GPSDB.pause 
hold  off 

clg 

cic 

%  For  comparison  purpose  we  plot  the  BPSK  signal  calculated 
%  without  interference  and  the  theoretical  one. 

axis([fc-5e6  fc+5e6  -90  0]); 
plot(f.PSSl(l;16384))... 
tiUeCBPSK  SIGNAL  COMPARISON  )... 
xlabeK'FREQUENCY,  Hz’),.. 
ylabcK’POWER.  dB'),.. 

text(0.13,0.85.' . THEORETICAL  (dotted  line)’,.. 

text(0.13,0.80,' _  CALCULATED  (solid  line)',.. 

pause 
hold  on 

plot(f.Gl(l:16384))... 

pause 

hold  off 

clg 

clc 

save  temp.FFT  SD  f  FIL  Pssl  S  fc 
end 


BLOCK  Signal  Filtering  on  the  Detector  Array  (S_FIL.m) 


echo  on 
load  temp_FFT 

clc 

%  Now  is  time  to  pass  the  BPSK  signal  through  the  detector  array 

SF=SD.*FIL;  %  Filtering  the  Signal 

Psf=SF.*conj(SF); 

Psfl=Psf/max(Psf);  %  PSD  Filtered  Signal 

clc 

%  First  of  all  we  compare  the  signal  spectrum  with  the  filter  due 
%  to  the  detector  array 

axis([fc-6e6  fc+6e6  0  1.2]); 
plot(f,Pssl(l;  16384), f.FILd;  16384).’;')... 
titleCFILTER  SIGNAL  &  BPSK  SPECTRUM  ')... 
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xlabeK'FREQUENCY,  Hz'),.. 
ylabelC AMPLITUDE,  volts'),., 
meta  GSFIL, pause 


cig 
c  Ic 

%  Now  we  can  plot  the  signal  filtered  compared  to  the  original  one 

axis([fc-6e6  fc+6e6  0  1.2]); 
pIot(f.Psfl(l:16384).f,Pssl(l;16384).':')... 
lilleCBPSK  SPECTRUM  COMPARISON’)... 
xlabeK’FREQUENCY.  Hz')... 
ylabeK’AMPLITUDE.  volts’)... 
pause 

clg 

cic 

save  temp_SF  SF  f 
delete  temp_FFT.mat 
end 


BLOCK  IFFT  Signal  Filtered  (IFFT_B,m) 


echo  on 
load  temp_SF 
load  temp_S 

c  Ic 

%  Now  we  return  to  the  time  domain  taking  the  inverse 
%  Fourier  transform  of  the  filtered  signal 

rs=ifft(SF);  %  Inverse  Fourier  transform 

t=t';  %  Arrange  time  axis 


c  Ic 

%  We  can  plot  the  Recovered  Signal  in  the  time  domain 

axis([0  0.00002  -Ac-1  Ac+lj) 
plot(t.rs),title(’RECOVERED  BPSK  SIGNAL’)... 
xlabeK’TIME.  sec’),.. 
ylabeK'AMPLITUDE,  volts’)... 
pause 

clg 

clc 

%  Let’s  compare  with  the  original  signal  plotting  it 
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subplot(211),ploi(t,s);  lille(’ORIGINAL  BPSK'),.. 
xIabeK'TIME.  sec').ylabel(’AMPLITUDE')... 
subplot(2l2),ploi(t.rs);tille('RECOVERED  BPSK  SIGNAL'),.. 
xIabeK’TIME,  sec’),ylabel('AMPLITUDE'),.. 
pause 

clg 

cic 

save  temp_R  rs  t  s  m  Ac  fc 
delete  iemp_SF.mat 
delete  temp_S.mat 
end 


BLOCK  Coherent  Detector  (COD_B.m) 


echo  on 
load  temp_R 
load  temp_SR 
cIc 

%  Now  is  time  to  make  a  coherent  detection.  The  optimum 
%  detection  consist  of  product  detection  with  an  integrator 
%  that  is  turned  on  at  the  beginning  of  each  bit  period 
%  and  integrates  to  the  end  of  the  period. 


DET=rs,*cos(2'*pi*fc'^t); 

% 

Coherent  detection 

B=ones(  1  ,n); 

% 

Base  vector 

for  u=0:m;n-m; 

% 

Loop 

sa=cumsum(DET(  1  ,(u+l  :u+m))); 

% 

Integrator 

B(  1  ,(u+ 1  ;u+m))=B(  1  ,(u+ 1  ;u+m)).*sa; 

sa=0; 

% 

Dump 

end 

% 

End  loop 

P=-l; 

H=ones(  1  ,n); 

for  v=0:m:n-m; 

% 

Loop 

dig=B(I,v+m/2); 

% 

Reading  the  middle 
point 

if  dig>0 

H(  1  ,(v+ 1  ;v+m))=H(  1  ,(v+ 1  ;v+m)).*  1 ; 

% 

Mapped  into  a  digital 
value 

elseif  dig<0 

H(  1  ,(v+ 1  ;v+m))=H(  1  ,(v+l  ;v+m)).*P; 

% 

Mapped  into  a  digital 
value 

end 

% 

End  conditional 
statement 

end 

% 

End  Loop 
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E=A.*H; 


Checking  for  error 


c  Ic 

9c  We  look  ai  ihe  output  of  the  integrator 

axis([0  a  -m-1  m+1]); 
plotd.B),.. 

titleCINTEGRATOR  OUTPUT'),.. 
xlabcU'TIME,  sec'),.. 
ylabeU'AMPLITUDE,  volts'),.. 
pause 

clg 

clc 

%  Comparison  of  the  Random  Binary  Data  and  the  data 
%  recovered 

axis([0  a  -Ac-2  Ac+2]); 

subplot(21  l),plot(t,A);  lille('ORIGINAL  BINARY  DATA'),.. 
xlabelCTIME,  sec'),ylabel('AMPLITUDE'),.. 
subplot(212),plot(t.H);title('RECOVERED  BINARY  DATA')... 
xlabeK'TIME,  sec’).ylabei(’AMPLITUDE'),.. 
pause 

clg 

clc 

9c  In  order  to  determine  if  an  error  is  present  we  compare 
9c  the  data.  If  we  get  all  I,  there  is  "no  error";  othcnvi.se  a  -I 
9c  represents  an  "error". 

axis([0  a  -2  2)); 

plot(t,E),title('ERROR  DETECTION'),.. 
text(0.5,0.75,'NO  ERROR','sc')... 
text(0.5,0.25,'ERROR','sc’),.. 
xlabeK'TIME,  sec'),.. 
ylabelCAMPLITUDE,  volts')... 
pause 

clg 

clc 

delete  temp_S.mat 
delete  temp_R.mat 
end 
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APPENDIX  C  -  INFORMATION  ABOUT  NOTCHES 


In  order  to  help  the  user  in  the  selection  of  the  interference 
frequency  and  evaluate  the  effect  in  a  particular  notch,  here  is  a 
table  showing  the  position  of  each  notch  depending  how  wide  it  is  in 
the  frequency  domain. 


Notch 

Notch 

Notch 

10  KHz 

Notch 

20  KHz 

# 

Center 

Start  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65895000 

65905000 

65890000 

65910000 

2 

66650000 

66645000 

66655000 

66640000 

66660000 

3 

67400000 

67395000 

67405000 

67390000 

67410000 

4 

68150000 

68145000 

68155000 

68140000 

68160000 

5 

68900000 

68895000 

68905000 

68890000 

68910000 

6 

69650000 

69645000 

69655000 

69640000 

69660000 

7 

70400000 

70395000 

70405000 

70390000 

70410000 

8 

71 150000 

71145000 

71155000 

71140000 

71 160000 

9 

71900000 

71895000 

71905000 

71890000 

71910000 

10 

72650000 

72645000 

72655000 

72640000 

72660000 

1 1 

73400000 

73395000 

73405000 

73390000 

73410000 

12 

74150000 

74145000 

74155000 

74140000 

74160000 

Notch 

Notch 

Notch 

30  KHz 

Notch 

40  KHz 

# 

Center 

Start  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65885000 

65915000 

65880000 

65920000 

2 

66650000 

66635000 

66665000 

66630000 

66670000 

3 

67400000 

67385000 

67415000 

67380000 

67420000 

4 

68150000 

68135000 

68165000 

68130000 

68170000 

5 

68900000 

68885000 

68915000 

68880000 

68920000 

6 

69650000 

69635000 

69665000 

69630000 

69670000 

7 

70400000 

70385000 

70415000 

70380000 

70420000 

8 

71150000 

71135000 

71165000 

71 130000 

71170000 

9 

71900000 

71885000 

71915000 

71880000 

71920000 

10 

72650000 

72635000 

72665000 

72630000 

72670000 

1  1 

73400000 

73385000 

73415000 

73380000 

73420000 

12 

74150000 

74135000 

74165000 

74130000 

74170000 

Notch 

Notch 

Notch 

50  KHz 

Notch 

60  KHz 

# 

Center 

Start  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65875000 

65925000 

65870000 

65930000 

2 

66650000 

66625000 

66675000 

66620000 

66680000 

3 

67400000 

67375000 

67425000 

67370000 

67430000 

4 

68150000 

68125000 

68175000 

68120000 

68180000 

5 

68900000 

68875000 

68925000 

68870000 

68930000 

6 

69650000 

69625000 

69675000 

69620000 

69680000 

7 

70400000 

70375000 

70425000 

70370000 

70430000 

8 

71150000 

71125000 

71175000 

71120000 

71180000 

9 

71900000 

71875000 

71925000 

71870000 

71930000 

10 

72650000 

72625000 

72675000 

72620000 

72680000 

1 1 

73400000 

73375000 

73425000 

73370000 

73430000 

12 

74150000 

74125000 

74175000 

74120000 

74180000 

Notch 

Notch 

Notch 

70  KHz 

Notch 

80  KHz 

# 

Center 

Start  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65865000 

65935000 

65860000 

65940000 

2 

66650000 

66615000 

66685000 

66610000 

66690000 

3 

67400000 

67365000 

67435000 

67360000 

67440000 

4 

68150000 

68115000 

68185000 

68110000 

68190000 

5 

68900000 

68865000 

68935000 

68860000 

68940000 

6 

69650000 

69615000 

69685000 

69610000 

69690000 

7 

70400000 

70365000 

70435000 

70360000 

70440000 

8 

71150000 

71115000 

71185000 

71110000 

71190000 

9 

71900000 

71865000 

71935000 

71860000 

71940000 

10 

72650000 

72615000 

72685000 

72610000 

72690000 

1 1 

73400000 

73365000 

73435000 

73360000 

73440000 

12 

74150000 

74115000 

74185000 

74110000 

74190000 

Notch 

Notch 

Notch 

90  KHz 

Notch 

100  KHz 

# 

Center 

Start  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65855000 

65945000 

65850000 

65950000 

2 

66650000 

66605000 

66695000 

66600000 

66700000 

3 

67400000 

67355000 

67445000 

67350000 

67450000 

4 

68150000 

68105000 

68195000 

68100000 

68200000 

0 

68900000 

68855000 

68945000 

68850000 

68950000 

6 

69650000 

69605000 

69695000 

69600000 

69700000 

7 

70400000 

70355000 

70445000 

70350000 

70450000 

8 

71150000 

71 105000 

71195000 

71100000 

71200000 

9 

71900000 

71855000 

71945000 

71850000 

71950000 

10 

72650000 

/ 2605000 

72695000 

72600000 

72700000 

1  1 

73400000 

73355000 

73445000 

73350000 

73450000 

12 

74150000 

74105000 

74195000 

74100000 

74200000 

7  5 


Notch 

Notch 

Notch  110  KHz 

Notch 

120  KHz 

# 

Center 

Stan  freq. 

End  freq. 

Start  freq. 

End  freq. 

1 

65900000 

65845000 

65955000 

65840000 

65960000 

2 

66650000 

66595000 

66705000 

66590000 

66710000 

3 

67400000 

67345000 

67455000 

67340000 

67460000 

4 

68150000 

68095000 

68205000 

68090000 

68210000 

5 

68900000 

68845000 

68955000 

68840000 

68960000 

6 

69650000 

69595000 

69705000 

69590000 

69710000 

7 

70400000 

70345000 

70455000 

70340000 

70460000 

8 

71150000 

71095000 

71205000 

71090000 

71210000 

9 

71900000 

71845000 

71955000 

71840000 

71960000 

10 

72650000 

72595000 

72705000 

72590000 

72710000 

1 1 

73400000 

73345000 

73455000 

73340000 

73460000 

12 

74150000 

74095000 

74205000 

74090000 

74210000 

To  put 

the  correct 

value  in 

the  notch  ! 

filter  that  represent  a  given 

attenuation  in  dB,  a  table  with  different  values  of  attenuation  from  0 
to  -45  dB  is  presented. 


ATTENUATION  TABLE 


Attenuation 

Number 

Attenuation 

Number 

(dB) _ 

_ (dB) 

0 

1.00000 

-11 

0.07943 

- 1 

0.79432 

-12 

0.06309 

-2 

0.63096 

-13 

0.05119 

-3 

0.50119 

-14 

0.03981 

-4 

0.39811 

-15 

0.03162 

-5 

0.31623 

-20 

0.01000 

-6 

0.25119 

-25 

0.00316 

-7 

0.19953 

-30 

0.00100 

-8 

0.15849 

-35 

0.00032 

-9 

0.12589 

-40 

0.00010 

-10 

0.10000 

-45 

0.00003 

APPENDIX  D  -  DEFINITIONS 


A  group  of  useful  definitions  related  to  this  work  are  presented. 

Bragg  cell  :  The  heart  of  all  acousto-optic  systems.  A  thin  slab 
of  a  transparent  crystal  (such  as  tellurium  oxide,  lithium  niobite, 
or  gallium  phosphide)  with  one  or  more  piezoelectric  transducers 
attached.  When  a  microwave  signal  excites  the  transducer,  it 
creates  an  ultrasound  wave  into  the  crystal.  This  wave  changes 
the  optical  index  of  refraction  at  the  microwave  frequency  [Ref. 
9]. 

Excisor-detector:  An  array  of  photodetectors  which  is  inserted 
into  the  frequency  plane  of  an  optical  spectrum  analyzer.  The 
small  spaces  between  detectors  are  used  as  spatial  light 
modulators,  that  produce  notches  in  the  frequency  response. 

Model  :  Mathematical  representation  of  a  system  or  part  of  a 
system  and/or  some  portion  of  its  environment.  The  model  is  a 
set  of  equations  which  may  be  solved  manually  or  by  means  of  a 
computer  to  determine  the  behavior  of  the  system  under  a 
specific  set  of  conditions  or  during  a  series  of  sets  of  conditions 
Models  permit  easy  modification  of  the  system  parameters  to 
determine  the  effects  of  such  changes  on  system  performance 
[Ref.  10]. 
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Simulation  :  Simulations  include  techniques  for  representing  a 
system  or  part  of  a  system  and/or  some  portion  of  its 
environment  (including  human  interfaces  with  hardware) 
through  the  use  of  sustitutes  for  all  or  part  of  the  real  hardware, 
personnel;  and  environment.  Thus  simulation  covers  a  more 
extensive  field  than  modeling,  although  most  models  are  also 
simulations  [Ref.  10]. 
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